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PREFACE 


This  paper  was  written  expressly  for  Professor  Henry  T. 
Nagamatsu,  Professor  of  Aeronautical  Engineering,  Rensselaer 
Polytechnic  Institute,  Troy,  New  York.  The  length 
limitation  of  the  paper  prohibited  a  detailed  discussion  on 
the  elementary  principles  of  shock  tube  operation  and  other 
gas  dynamic  phenomena.  For  this  reason,  readers  unfamiliar 
with  these  concepts  should  first  begin  with  a  thorough 
review  of  the  principles  of  unsteady,  one  dimensional  gas 
dynamics.  An  excellent  text  for  this  is  Gas  Dynamics.  2nd 
ed.  by  James  E.  A.  John,  Allyn  and  Bacon,  Inc.  Newton,  MA, 
1984.  Also,  readers  unfamiliar  with  the  elements  of 
computational  fluid  dynamics  should  begin  with  a  review  of 
these  elements.  An  excellent  text  for  this  is  Reference  1. 

I  am  extremely  grateful  to  Dr.  Gary  Carofano  of  the  US 
Army’s  Benet  Weapons  Laboratory  for  the  work  he  did  in 
applying  Harten’s  scheme  to  an  easily  workable  computer 
code.  In  addition,  his  many  informative  and  stimulating 
discussions  led  to  an  easier  understanding  of  the 
complications  of  CFD  limitations,  the  thermodynamics  of  the 
problem,  and  an  overall  better  understanding  of  these 
phenomena.  His  patience  and  depth  of  understanding  are 
greatly  appreciated. 
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A  iiL.'  TRACT 


A  auraer  i-  a study  was  inaue  cl  the  t.p  time  in  a 
hvpf-rionic  .ho-  k  tunnei  tor  Ma  -  h  10  flow.  liar  ten*  ■=»  totcil 
variation  d  i  mi  n  ish  i  ng  ‘TVUi  heme  was  used  to  simulate  t  he 
entire  shock  tunnel  operation  from  diaphragm,  rupture  to 
disruption  oi  the  1  low  in  the  nozzle.  The  -.ode  was  able  to 
.r.odel  the  main  features  of  unsteady  flow  in  a  •-,ho«~k  tube  to 
include  the  shock  wave  and  contact  surface  Intersections. 
Results  were  obtained  that  predict  the  steady  flow  duration 
time  for  stagnation  temper iture =  of  1200’F,  2700'  F,  and 
o'uOO  ’F.  Both  air  and  helium  were  used  as  driver  gases  with 


Ci ui  ijl  e r  One 


NTROUL  .7  ION 


OBJ  ECT  I  VE  OF  7  HE  STUDY 


This  work  was  undertaken  to  accomplish  two  goals:  to 

conduct  a  numerical  analysis  of  the  test  time  available  in 
the  Rensselaer  Polytechnic  Institute  (RPI )  hypersonic  shock 
tunnel  arid  determine  the  conditions  yielding  the  maximum 
r>ossi  ble  flow  duration  time. 


HISTORICAL  BACKGROUND 


In  the  1st  e  IGGOs  and  early  1960s  much  aeronautical 
research  was  being  devoted  1o  tne  high  Mach  number  regime 
known  as  hypersonic  (27:1).  The  research  centered  primarily 
oil  experimental  <3,  17-.2G,  32, 34  -:_,9,  4  1 , 42)  and  analytical 
(10,26,20,30)  work  in  support  of  such  programs  as  Mercury, 
Gemini,  Apollo,  and  X-1‘5  for  the  National  Aeronautics  and 
Space  Administration  (NASA/  arid  research  on  ballistic 
missile  reentry  vehicles  for  the  Air  Force  and  Navy  (44:2). 
In  the  early  1970s,  the  X- 1  b  and  Apollo  projects  were 
completed,  America’s  Super  conic  Transport  oSST)  was  killed 
by  Congress,  arid  the  bulk  of  the  design  work  for  the  space 
shuttle  was  completed.  Few  resell ch  dollars  were  available 
for  hypersonic  aerodynamics  causing  the  effort  in  this  area 
to  stop  :44:3).  Virtually  ail  of  the  hypersonic  experimental 
research  facilities  in  tin-  United  States  fell  into  a  state 
of  disuse  or  were  dismantled  and  sold  for  scrap  '27:1). 

Until  1984  ii>  p-*i  sorri  c  research  In  the  United  States  was 
nearly  none::  i  stent .  List.,  of  some  of  the  hvpei  sonic 
research  to  da  1  e  ms  given  In  References  27  and  33. 


Today  the  i  *.•  rearch  In  hvpei  ionics  is  mainly  being 
conducted  in  the  growing  field  of  computational  fluid 
dynamics  ;i'.FL/  .  This  was  made  possible  by  the  tremendous 
increases  in  the  memory  and  •:  omputat  i  ona  1  power  of  modern 
supercomputers.  I +.  is  now  pos.-.ible  to  obtain  the  complete 
three-dimensional  inviscid  flow  field  around  a  hypersonic 
vehicle.  CFD  has  the  advantages  of  being  aide  to  yield 
detailed  pictures  of  the  1  low  iields  uncontaminated  by 
intrusive  measuring  device',  which  disturb  the  flow.  The 
ia.ge  expense  of  exper  i  men  ‘  a  i  work  makes  CFD  ruore  attract  ive 


n.o t  Hi  i  vt-  r  i  .  i  e -s  . 


The  Cr  •  1 1  input  ei  •/ode--,  h;  iwvvr1 1  are 
only  approximations  oi  the  phyc.i  ■-  e  ot  the  flow  and  i  equiiw 
?•  :  rap  l  i f y  Lug  as:iUH|)t  lone  to  keen  memory  requirements  and 
computational  times  reasonable  even  with  the  mus-  powerful 
super  compute  r  s .  There  are  some  expo  imental  fa'  ilifies 
still  being  used  however.  In  order  to  initially  validate 
the  C?t>  methods,  experimental  results  must  be  ir-d .  This  i.; 
t  h*.-  only  way  to  one.-,  k  the  validity  of  t.he  computational 
re  ..nits  when  dealing  with  the  real  gas  effects  ot  hyper  son!', 
flows. 


SIGNIFICANCE  U- 


PROBLEM 


In  ,  the  interest  in  hyper  .onlcs  was  prompted  by  tin 

need  for  niore  knowledge  to  support  new  vein  ole  concepts. 


of  these  project: 


«  i  '  The  Aero  i  i  sted  Orbital  Transfer-  Vehicle 
wVliyj.,  to  be  used  to  tr.5iis.ter  payloads  between 
g  e  os  y  r  i  hronuus  orbit  to  »  ow  earth  orbit  <  LEO  >  to  a 
rendezvous  with  the  shuttle-  during  the 
aerodynamic  braking  stage,  the  vehicle  will  travel 
hyper  son  ica  1  iy  (Mach  30  >  through  the  earth’s 
a  t  mc'-.phet  • -  at.  a'titudes  ot  ."-i'VgOOO  tn  fsOO.t'OG  feet. 


-  '  The  National  Aerospace  Plane  <■  NASP >  .  This 
Air  Force  vehicle  will  take  oil  and  land 
not  r.-.ontaily  from  convent  1  otia  i  airports,  be  able 
to  travel  at  Mach  12,  go  into  orbit,  deorbit  and 
land.  The  civilian  version  of  this  vehicle  will 
crave)  at  about  Mach  7  and  about  110,000  feet.  This 
is  the  so-called  Orient  Express  Hyper  sonic 
Tr  aiispor  t  ’  HoT  >  . 


•id)  A  sec end -gone rat  ton  space  shuttle.  As 
technology  arid  knowledge  increase,  a  second-genera  t  i  on 
shrift  le  will  n.-.  J  to  be  built  to  fulfill  t  tie  need 
for  a  follow  on  space  lawn,  h  vehicle. 


■  4  ■  Low-hypersonic  missl  lc  anti-tank  penet  txiturs. 

T.ie  Army  will  use  hypersonic  <  Mach  5-7),  tin-stabilised, 
armor  -  piercing  projectile:-,  to  penetrate  advanced  armor? 
on  future  tank-,  and  armored  vehicles. 


1  ,  w 


1  i  i*r  a. 

i  e  i  d  -j 


i  1  I  u  -  ■ 1  r  a  i-  e  ►  h»*  c  h.  •  i  letige  i  n  studying  hypersonic 
,  consider  a  typical  i  i  I’  M  wi  t  h  a  'moo  mi  le 
h  j  v-.'ihi'  1  e  reenter-,  the  a t mosph*-" t  e  ,it  Mach  2u. 

.-I  t.he  air  t  e  rape  i  a  t  1 1  r  e  i ..  >h  i  nd  I  lie  n  f  i  r  m.  i  1  .lit  >c  k 
:  bo  I  v  .  in  ‘  a  .  t.  ign  a  oOnuk'  >  !<■,  it'  pi 


.’.l  .  o  n 


<17:  34  1  >■  A  •  f  his  t  e  ru  |  >e  r  a  i  hi  e  lit  no  Innj/M  b.-have^ 

I  <J.-»cj  J  gas  ml  Lure  ut  <  i  i-t  1  <>  in  i  •  oxygon  ami  n  i  t  f  ugen .  I  nst  »-.<■: , 
j!  wi  i  ]  .tKidl.)  a  in  1  ion  i z.e.  The  resulting  "real  gas" 

will  i.on&isl.  Cif  monatomic  ami  ionic  oxygen,  nitrogen,  oxide" 
of  nitrogen,  and  free  electrons.  As.  this-,  real  gas  mixture 
travels,  bark  along  the  vehicle  away  from  the  stagnation 
point,  it.  expands  and  coole.  There  may  be  .oinc  chemical 
recombiriatlc.il  if  the  press-ire  is  high  enough.  If  not,  the 
original  mi:- line  will  be  "1  ru<:t:i"  at  its  original 
composition.  These  non-equilibr  iuni  effects  cause  a  large 
departure  from  conventional  aerodynamic  and  heat  transfer 
theories.  These  effects  must,  be  fully  understood  and 
accounted  for  in  the  design  stage  for  any  hypersonic 
vehicle.  The  only  way  to  study  these  effects  is  to 
duplicate  the  hypersonic  flow  conditions  in  the  laboratory. 

GROUND  TEST  FACILITIES 

The  experimental  fa.ilities  used  today  for  hypersonic 
research  are  of  three  basic  types:  continuous,  blowdown,  and 
impulse.  Table  1  gives  the  characteristics  of  each  type. 
Conventional  experimental  tools  for  getting  aeronautical 
information,  such  as  supersonic  wind  tunnels,  are  inadequate 
for  invest  igai.  ing  the  real  gas  effects  of  the  flow  around 
vehicles  at  hypersonic  speeds  because  of  the  reservoir 
temperature  limitations.  An  impulse  facility  (shock  tubes, 
shock  tunnels,  and  gun  tunnels.)  is  the  only  facility  that 
can  exactly  duplicate  hypersonic  flow  conditions.  For  this 
reason  some  researcher s  continue  to  use  shock  tunnels  for 
experimental  hypersonic  research.  (For  an  excellent 
description  on  the  operation  of  hypersonic  shock  tunnels  see 
Reference  30.  >  Figure  1  gives  a  good  description  of  the  test 
times  and  maximum  stagnation  temperatures  attainable  in  the 
various  ground  test,  facilities.  Since  the  flow  duration  in 
a  shook  tunnel  is  between  1  and  10  milliseconds  (ms),  it  is 
imperative  that  the  conditions  be  optimized  for  each  run  to 
obtain  the  maximum  test  time  possible.  Note  that  a  test 
time  of  one  millisecond  is  not  as  limiting  as  it  may  sound. 
The  response  times  of  pressure  transducers  and  thin  film 
heat  gauges  at e  on  the  order  of  a  microsecond.  In  order  to 
aid  the  experimenter  in  deter min ing  the  best  possible  test 
times,  in  a  hypersonic  shock  tunnel,  CFD  can  be  used  to  model 
Uie  starting  process  in  th~  sho.k  tube  and  the  flow 
processes  in  the  hypersonic  nuzzle  and  test  section.  This 
is  easier  and  can  yield  more  accurate  results  than  running 
several  sets  of  conditions  ir  the  shock  tunnel  to  determine 
the  best  test  times. 


‘H 

a 

?  *i 


1 

Continuous  Flow 

Blow  Pawn 

I mpu lse 

Mach  No. 

1.5  to  10 

i.4  -  14 

7-86 

Run  Time 

cant i nuous 

seconds  to 
inf  nutes 

1  -  100  ms 

Test  Time 

long 

med  i  urn 

short 

Repeatabl 1 ity 

excel  lent 

medl  urn 

1  ow 

Cost 

very  high 

high 

1  ow 

Stagnat ion 
Temperature  T. 

,  low 

i  ow 

high 

Set-up.  Time 

very  long 

shor  t 

short 

Reynolds  No. 

low 

medium 

high 

Instrumentation  standard 

standard 

spec ial 

Productivity 

very  high 

medium 

low 

TABLE  1.  Comparison  of  test  facilities  (45:7,8,10). 


The  test  time  (flow  duration  time)  in  a  hypersonic  shock 
tunnel  is  limited  by  the  starting  process  in  the  nozzle,  and 
the  passage  of  either  the  reflected  sho<"  k  from  the 
shock- contact  surface  intersection  or  the  reflected  head  of 
the  expansion  wave  from  the  shock  tunnel  driver  section 
(30:--;  40:--).  The  aim  of  the  experimenter  is  to  achieve 
the  maximum  possible  flow  duration  time  with  a  high  enough 
stagnation  temperature  (high  incident  shock  Mach  number)  to 
duplicate  hypersonic  reentry.  A  sample  wave  diagram  is 
shown  in  Figure  2  to  show  the  significance  of  the  test  time. 


The  shock  tunne 
rupture  of  the  diap 
driven  tube  until  i 
shoulders  and  trave 
also  travels  into  t 
nozzle,  this  starti 
and  a  rearward  faci 
high  velocity  gas  f 
once  this  swallowed 


1  produces  a  strong  shock  wave  by  the 
hragra.  This  shock  wave  travels  down  the 
t  partially  reflects  off  the  nozzle 
Is  toward  the  driver  section.  The  shock 
he  nozzle  starting  the  flow.  In  the 
tig  shock  is  followed  by  a  contact  surface 
ng  shock  that  is  "swallowed”  by  the 
low  in  the  nozzle.  Testing  can  commence 
shock  clears  the  test  section  <40:2). 


Figure  1.  Stagnation  temperature  and  flow  duration  time  for  test  facilities  (45;  Figure  6). 


Figure  2.  Typical  wave  diagram  for  a  hypersonic  shock 
tunne 1 . 


Iii  the  roost  <-  onunou  case  u*:  <  ui  ring  In  d  shock  tunnel,  the 
reflected  shock  eventually  intersects  the  approaching 
contact  surface.  For  most  cases  this  intersection  produces 
two  shocks,  one  that  is  transmitted  through  the  contact 
surface  toward  the  driver,  and  a  reflected  shock  that 
travels  back  toward  the  nozzle.  The  reflected  shock  then 
partially  reflects  off  the  nozzle  shoulders  and  intersects 
the  contact  surface  again.  It  also  continues  into  the 
nozzle  where  it  ends  testing  when  it  disturbs  the  flow 
conditions  in  the  nozzle  (45:4).  Two  other  cases  may  occur. 
One  is  the  unsteady  expansion  head  may  reflect  off  the 
driver  end  wall,  overtake  the  contact  surface  and  travel 
into  the  nozzle.  This  ends  the  test  time  before  any  shock 
wave-contact  surface  intersection  can  take  place  This  is 
more  common  with  very  short  drivers,  gases  with  high 
acoustic  velocity  in  the  driver  (such  as  helium  or 
hydrogen),  and  very  low  shock  Mach  numbers.  The  other  case 
occurs  if  a  gas  with  a  high  acoustic  speed  is  used  in  the 
driver  (helium,  hydrogen  or  hydrogen/oxygen  combustion 
drivers).  This  will  produce  a  "soft"  contact  surface 
(4:446-447).  When  the  reflected  shock  intersects  this  soft 
contact  surface  it  produces  a  transmitted  shock  but  instead 
of  a  reflected  shock  a  reflected  expansion  is  produced.  The 
head  of  this  expansion  travels  into  the  nozzle  at  the  local 
fluid  velocity  plus  the  local  speed  of  sound. 
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Chapter  Two 


NUMERICAL  METHODS  FOR  HYPERSONIC  FLOW 


ASSUMPTIONS  AND  LIMITATIONS 

Two  major  factors  in  choosing  a  numerical  method  were 
the  computing  power  available  and  the  applicability  of  a 
particular  scheme.  The  computer  used  for  this  study  was  a 
Prime  850  minicomputer — a  relatively  slow  machine  with 
virtual  memory  architecture.  The  code  was  kept  simple  to 
preclude  excessive  computing  times.  This  was  the  major 
reason  for  rejecting  the  CFD  methods  which  account  for  fluid 
viscosity  and  real  gas  effects.  Computing  was  also  limited 
to  one-dimensional  where  passible.  The  only  place  in  the 
shock  tunnel  where  two-dimensional  computations  were  done 
was  in  the  nozzle  section  where  strong  two-dimensional 
effects  were  found.  A  look  at.  the  literature  shows  some 
numerical  work  centered  mainly  on  the  nozzle  starting 
process  but  not  on  the  entire  shock  tunnel  operation  <40: 2>. 
An  excellent  text  on  CFD  is  Reference  1.  This  gives  several 
numerical  methods  that  are  app] icable  for  hypersonic  flow 
and  this  study.  Reference  40  describes  previous  work  on 
the  nozzle  starting  process,  while  Reference  15  describes 
the  propagation  of  a  shock  through  the  nozzle.  The  computer 
code  used  in  Reference  40  was  a  good  first-order  scheme  but 
was  found  unsuitable  for  this  study  due  to  its  tendency  to 
smear  the  shocks  somewhat  and  the  contact  surfaces  greatly. 


DESCRIPTION  QF  NUMERICAL  METHOD  USED 

One  of  the  most,  critical  interactions  affecting  the  test 
times  in  hy per  conics  shock  tunnels  is  the  interaction  of  the 
reflected  shock  and  the  contact  surface.  If  the  contact 
surface  is  not  properly  treated,  inaccurate  estimates  of  the 
test  time  could  result.  Foi  thi-  reason  a  numerical  method 
which  has  good  performance  in  treating  contact  surfaces  was 
sought.  Such  a  method  was  developed  and  refined  by  Harten 
<10: — )  and  written  into  a  computer  code  for  this  study  bv 
Carafano  (28.--).  The  Harden  method  does  not  smear  the 
contact,  surface  like  the  first  order  methods  such  as 
Godunov,  Roe,  or  Lax.  It  also  does  not  have  the  instability 
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problem  near  severe  discontinuities  like  the  second-order 
methods  of  MacCorraack,  Lax-Vendrof f ,  and  others.  For  these 
reasons  the  Harten  method  describe-’  In  References  10  and  28 
was  chosen  for  this  study. 

The  Harten  method  is  a  second-order  accurate  numerical 
scheme  which  solves  the  unsteady  compressible  Euler 
equations,  the  continuity  equation,  and  the  energy  equation. 
The  code  as  written  treats  the  fluids  as  perfect  gases  with 
no  viscosity  (28:13-15).  The  code  takes  into  account  two 
separate  gases— one  for  the  driver  and  one  for  the  test  gas. 
Thus  the  code  can  handle  helium  driven  shock  tunnels  as  well 
as  combustion  driven  tunnels.  However,  combustion  driven 
shock  tunnels  should  be  handled  wii h  caution.  Since  the 
code  only  handles  perfect  gases,  the  maximum  shock  Mach 
number  must  be  kept  below  about  six  to  preclude  the 
occurrence  of  real  gas  effects  (dissociation  and  ionization; 
caused  by  the  very  high  stagnation  temperatures. 

The  cade  treats  the  shock  tube  portion  of  the 
calculations  as  strictly  one-dimensional.  The  hypersonic 
nozzle  can  be  treated  as  either  one-dimensional  or 
two-dimensional  axi-symmetr ic .  Graphics  packages  were 
written  that  produce:  two-dimensional  contour  plots  of  the 

flow  in  the  nozzle,  history  plots  of  thermodynamic 
properties,  and  distributions  >f  pressure  and  Mach  number  in 
tiie  shock  tube. 

The  RPI  shock  tunnel  modeled  in  this  study  is  shown  in 
Figure  3.  An  enlargement  of  the  nozzle  area  is  shown  in 
Figuj e  4.  The  transducer  locations  are  shown.  These 
allow  histories  to  be  made  of  the  state  variables.  Table  2 
shows  the  exact  transducer  locations  as  measured  from  the 
driver  end  wall.  The  Mach  number  shown  is  the  Mach  number 
at  that  location  with  fully  developed  steady  flow. 


Cell  Position 

X  Position 
t Inches) 

Radius 
<  I  nches) 

A/A* 

Mach  No. 

1900 

810.  3 

2. 0 

16. 

0.03 

1950 

820. 

2. 0 

16. 

0.03 

2000 

829.3 

2. 0 

16. 

0.03 

2050 

838.  3 

2. 0 

16. 

0.03 

2100 

847. 

2.0 

16. 

0.03 

2110 

848. 5 

2. 0 

16. 

0. 03 

2140 

850.9 

0.5 

1. 

1.0 

2163 

852.6 

0.847 

2. 87 

2.59 

2183 

855.  4 

1 . 6 

10.3 

3. 95 

2203 

860.  7 

3. 0 

36.  6 

5.5 

2223 

870. 7 

5. 7 

129.  4 

7.  4 

2247 

894. 7 

12.  1 

587. 3 

10. 

Table  2.  Shock  tunnel  transducer  data. 
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Chapter  Three 


DISCUSSION  AND  ANALYSIS 
OF  NUMERICAL  DATA 


MACH  2.3  INCIDENT  SHOCK 


Air  Driver 


The  first  set  of  results  are  those  for  an  air  driven 
shock  of  Mach  2.3.  The  stagnation  temperature  obtained  for 
this  Mach  number  is  1200<:,f,  Other  Mach  numbers  and 
stagnation  temperatures  are  shown  in  Table  3. 


Incident  Mach  No. 

- TolTl - 

2.3 

1200 

3.41 

2700 

5.4 

6900 

Table  3.  Stagnation  temperature  and  shock  Mach  number. 


The  static  pressure  history  at  the  exit  of  the  nozzle  is 
shown  in  Figure  5.  The  top  trace  corresponds  to  the  Mach  10 
flow  conditions  at  the  exit  of  the  nozzle  and  the  bottom 
trace  is  at  the  Mach  7.4  position  inside  the  nozzle.  The 
ordinate  is  the  pressure  and  the  abscissa  is  the  time  in  ms 
since  diaphragm  rupture.  The  three  numbers  in  parentheses  to 
the  left  of  each  curve  are  the  transducer  cell  x  location,  y 
location  and  the  full  scale  value  of  the  plot.  All 
thermodynamic  properties  are  nondimensional ized  with  respect 
to  the  dump  tank  values. 

At  the  exit  of  the  nozzle  the  starting  shock  wave  is 
seen  as  a  rapid  rise  in  pressure.  The  rearward  or 
"swallowed”  shock  is  seen  as  the  rapidly  falling  pressure. 
The  resulting  steady  low  pressure  is  the  test  pressure  at 
the  exit  of  the  nozzle  indicating  fully  developed  Mach  10 
flow  conditions.  These  conditions  start  at  about  23.75  ms 


and  at  about  31.25  ms  a  rise  is  visible  in  the  pressure 
trace.  This  is  the  shock  wave  that  has  reflected  off  the 
nozzle  shoulders,  then  off  the  contact  surface,  and  traveled 
through  the  nozzle  disturbing  the  flow  conditions  and  ending 
the  test  conditions.  The  total  test  time  for  this  case  is 
about  7.5  ms.  A  plot  of  density  and  probe  total  pressure 
are  shown  in  Figures  6  and  7  respectively.  Additional  flow 
details  can  be  seen  in  these  plots  as  well  as  the  very 
well-defined  steady  state  test  conditions. 

The  termination  of  the  test  time  for  this  case  was 
caused  by  the  arrival  of  a  shock  wave  which  originated  with 
the  incident  shock  from  the  initial  diaphragm  rupture.  This 


PRESSURE 

TIME 
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Static  pressure  history  Mach  10  < upper!  and  Mach 
7.4  flower).  Mach  2.3  incident  shock,  air  driver. 


Figure  5. 


:: 


51 


| 


I 

1 

k*1 


Kw  "'S-'TX.^K^X.^V  .1  ST*  W  ^  v.  v%;  wvvtl'VV^'V.nr.'V.'7.',7  V 


shock  traveled  down  the  tube,  ruptured  the  nozzle  diaphra 
reflected  off  the  nozzle  shoulders,  traveled  back  upstrea 
toward  the  driver  end  of  the  tube,  reflected  off  the  contact 
surface  and  then  traveled  through  the  nozzle  to  the  test 
section.  Each  time  the  shock  hits  the  shoulders  of  the 
nozzle  a  transmitted  wave  travels  down  the  nozzle  and  a 
reflected  wave  travels  toward  the  driver.  Also,  each  time  a 
shock  wave  encounters  the  contact  surface,  a  transmitted 
wave  travels  toward  the  driver  and  a  reflected  wave  travels 
toward  the  nozzle.  This  complicated  wave  pattern  is 
sketched  in  Figure  2.  The  first  reflected  wave  is  the  only 
one  of  interest  in  this  study  since  it  is  responsible  for 


Figure  6.  Density  history  Mach  10  (upper)  and  Mach  7.4 

(lower).  Mach  2.3  incident  shock,  air  driver. 


a  oq 


ending  the  test  conditions  in  the  nozzle. 


Another  look  at  the  shock-contact  surface  interactions 
can  be  done  with  a  pressure  and  Mach  number  versus  distance 
plot  as  shown  in  Figure  8.  This  is  a  discrete  plot  oi 
static  pressure  and  local  Mach  number  distribution  along  the 
tube.  Figure  8  shows  the  distribution  14  ms  after  diaphragm 
rupture.  The  incident  shock  wave  can  be  seen  at  location 
o35 .  At  location  494,  the  contact  surface  is  visible  only 
in  the  Mach  number  plot,  since  the  pressure  is  equal  across 
a  contact  surface.  The  unsteady  expansion  can  be  seen  from 
about  location  352  to  the  driver  end  wall. 


Figure  O  shows  the  same  plot  at  a  later  time,  21  ms. 

The  shock  has  Just  ruptured  the  nozzle  diaphragm  at  location 
847.  The  contact  surface  and  expansion  are  evident.  It  is 
clear  from  the  two  plots  that  the  contact  surface  is  moving 
away  from  the  unsteady  expansion  which  will  result  in  the 
shock  reflecting  off  the  nozzle  shoulders  and  encountering 
the  contact  surface  and  not  the  expansion.  Figure  10  shows 
this  about  to  happen.  Notice  that  the  shock  has  reflected 
off  the  nozzle  shoulder  and  is  traveling  toward  the  contact 
surface.  The  reflected  shock  is  evident  from  the  high 
pressure  region  at  the  right  side  of  the  plot.  In  Figure 
11  the  shock  has  intersected  the  contact  surface  and 
resulted  in  a  reflected  shock  which  traveled  into  the 
nozz le . 

In  order  to  maximize  test  conditions  in  a  fixed  geometry 
shock  tunnel,  the  reflected  shock  from  the  contact  surface 
interaction  must  be  eliminated.  The  subject  of  shock  wave 
contact  surface  interaction  is  treated  very  well  in  several 
texts  <5: — ;  6:417-423).  This  can  be  accomplished  by 
matching  the  acoustic  impedance  of  the  gases  on  either  side 
of  the  contact  surface.  In  the  RPI  tunnel  this  can  be  done 
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by  using  helium  as  the  driver  gas  and  operating  at  the 
required  Mach  number  to  achieve  this  ’’tailored”  contact 
surface.  Using  helium  to  create  a  tailored  contact  surface 
calls  for  an  incident  shock  Mach  number  of  3.41.  Before  the 
results  of  the  tailored  case  are  presented,  the  results  of 
using  helium  and  a  shock  of  2.3  are  presented. 


Helium  Driver 


The  pressure  and  Mach  number  versus  distance  plot  for 
the  helium  driven  Mach  2.3  case  are  shown  at  28.5  ms  in 
Figure  12.  At  this  time  the  shock  has  already  reflected  off 
the  nozzle  shoulders.  Notice  the  sloping  traces  near  the 
right  side  of  the  plot,  This  indicates  that  the  reflected 
head  of  the  expansion  has  passed  through  the  contact  surface 
(visible  at  location  670)  and  has  traveled  into  the  nozzle. 
In  this  case  the  contact  surface  is  "soft"  and  an  expansion 
will  result  from  the  intersection  of  the  shock  and  contact 
surface  (6:420).  The  test  time  will  end  when  the  expansion 
moves  into  the  nozzle. 


Figure  9. 


Mach  number  and  pressure  distribution,  21.5  ms. 
Mach  2.3  shock,  air  driver. 
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The  lotcil  pressure  history  tor  this  case  is  shown  in 
Figure  13.  Notice  at  the  exit  of  the  nozzle,  flow  conditions 
last  only  about  0.5  ras .  This  is  much  shorter  than  the 
results  obtained  with  the  air  driver.  This  is  due  to  the 
acoustic  velocity  of  the  helium  being  about  four  times 
faster  than  that  of  air.  This  causes  the  reflected  head  of 
the  expansion  wave  to  close  on  the  shock  wave,  thereby 
cutting  test  time.  Also  the  reflected  shock-contact  surface 
intersection  yields  a  reflected  expansion  whi‘~h  moves  into 
the  nozzle  further  reducing  test  time.  This  type  of 
shock-contact  surface  interaction  is  shown  in  Figure  14 fe > . 
Therefore,  tor  low  Mach  numbers,  air  is  a  better  choice  for 
a  driver  gas  with  this  tube  layout. 


MACH  3.41  I NC I  DENT  SHOCK 

Using  room  temperature  gases,  a  helium  driver  and  air 
test  gas  can  produce  a  tailored  contact  surface  if  a  Mach 
3.41  incident  shock  is  used  (5:353).  With  this  setup  a 
reflected  shock  will  not  occur  when  a  shock  encounters  the 
contact,  surface.  These  flow  conditions  were  run  and  can  be 


moch  no. 
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seen  in  the  next  set  of  figures.  The  distribution  plot  is 
shown  for  this  case  in  Figure  15.  The  shock  has  already 
reflected  off  the  nozzle  shoulders,  has  passed  through  the 
contact  surface  without  reflecting  and  is  traveling  toward 
the  driver  end  wall.  The  sloping  pressure  curve  to  the 
right  of  the  shock  show  the  expansion  wave  already  passing 
through  the  shock  and  traveling  toward  the  nozzle.  Figure 
16  shows  the  probe  total  pressure  at  the  exit  of  the  nozzle. 
Notice  the  arrival  of  the  starting  shock  by  the  rapid  rise 
in  the  pressure  in  the  test  section,  the  "swallowed”  or 
rearward  facing  shock  by  the  sharp  drop  in  pressure  to  the 
steady  state  test  conditions.  The  test  conditions  terminate 
with  the  arrival  of  the  reflected  expansion  which  is  clearly 
seen  as  the  unsteady  decrease  in  pressure.  This  case 
resulted  in  4  ms  test  time.  This  would  not  have  been 
possible  with  air  as  the  driver  gas. 


MACH  5.4  INCIDENT  SHOCK 

For  the  last  case  of  this  study,  a  Mach  5.4  incident 
shock  corresponding  to  6900<'F  stagnation  temperature  was 
used.  At  these  conditions,  the  contact  surface  will  no 
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Figure  11.  Mach  number  and  pressure  distribution,  28  ms. 
Mach  2.3  shock,  air  driver. 


i 

I 


2 


MACH  NO. 


longer  be  tailored  and  a  "hard"  reflection  will  occur  such 
as  that  seen  in  the  air  driven  case  for  Mach  2.3.  Figure  17 
shows  the  pressure  at  the  last  two  transducer  locations.  At 
the  exit  of  the  nozzle,  the  test  conditions  last  for  about 
1.5  ms.  These  conditions  are  terminated  by  the  arrival  oi 
the  first  reflected  shock  from  the  contact  surface.  The 
Mach  5.41  shock  conditions  are  the  closest  to  the  Mach  10, 
150,000  ft.  cruise  conditions  of  about  7000"F  stagnation 
temperature . 
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Figure  12.  Mach  number  and  pressure  distribution,  28  ms. 
Mach  2.3  shock,  helium  driver. 
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Frobe  total  pressure  history,  Mach  10  (upper) 
and  Mach  7.4  (lower).  Mach  2.3  incident  shock, 
helium  driver . 
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Figure  14.  Wave  diagram  for  shock  interactions:  (a)  shock 
reflection  from  a  rigid  wall  or  closed  end;  <b)  shock 
incident  on  open  end  at  constant  pressure;  <c>  intersection 
of  shocks  of  opposite  family;  (d)  intersection  of  shocks  of 
common  family;  (e)  shock  incident  on  contact  surface;  and 
<f)  shock  incident  on  contact  surface  (6: Figure  8.35). 
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Figure  17.  Static  pressure  history  Mach  10  (upper)  and  Mach 
7.4  (lower).  Mach  5.4  incident  shock. 
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Chapter  Four 


CONCLUSIONS 


The  numerical  scheme,  Marten's  Total  Vai  iational 
Diminishing  (TVD>,  adapted  for  this  problem  by  Carofano 
<28:-->,  accurately  modeled  the  unsteady,  compressible  flow 
of  a  shock  tunnel.  The  TVD  scheme  also  accurately  modeled 
the  nozzle  starting  process  to  include  the  transversal  by 
the  starting  shock,  the  formation  and  passage  of  the  contact 
surface  and  the  rearward  facing  shock.  The  TVD  scheme  also 
accurately  modeled  the  shock  wave/contact  surface 
intersection  problem. 


Accurate  estimates  of  the  total  available  test  time  for 
the  RPI  hypersonic  shock  tunnel  were  obtained.  The  maximum 
test  time  occurred  for  the  Mach  2.3  incident  shock  with  air 
as  the  driver  gas.  This  time  was  found  to  be  7.5  ms.  For 
higher  stagnation  temperatures,  and  hence  higher  incident 
shock  Mach  numbers,  helium  is  needed  as  the  driver  gas. 

Under  these  conditions,  the  maximum  test  time  was  found  to 
occur  at  the  tailored  contact  surface  condition  or  Mach 
3.41.  The  test  time  at  this  condition  was  4  ms.  The  maximum 
test  time  occurred  with  air  as  the  driver  gas  and  low  shock 
Mach  numbers. 


At  the  highest  Mach  number  studied,  Mach  54,  the  test 
time  was  found  to  be  1.5  ms.  This  is  certainly  long  enough 
to  obtain  valuable  data  with  the  fast  response  time  of 
modern  instrumentation  and  signal  processing  techniques. 


Mach  No. 

Test  Time  (ms) 

Driver  Gas 

2! .  3 

7.5 

Air 

2. 3 

0.5 

Helium 

3.41 

4.0 

Hel ium 

5.  4 

1 . 5 

He  1 i um 

Table  4.  Comparison  of  Test  Times. 


Bl 


RECOMMENDATIONS 


Viscous  effects  need  to  be  studied.  The  test  times  can 
be  affected  by  the  formation  of  boundary  layers  and  other 
viscous  effects  and  hence  need  to  be  modeled. 

Adaptive  grid  techniques  need  to  be  employed  to  keep 
computation  times  down  and  to  obtain  higher  resolution.  The 
present  study  with  a  fixed  grid  and  mostly  one-dimensional 
computations  took  an  enormous  amount  of  computer  time  to 
solve  the  equations  within  acceptable  limits.  On  the  order 
of  30-45  CPU  hours  were  used  per  set  of  conditions. 

Lastly,  and  most  important,  experiments  need  to  be 
run  which  duplicate  the  numerical  simulations.  This  will 
allow  a  good  comparison  to  be  made  between  experimental  and 
numerical  results  to  determine  the  effectiveness  of  the 
present  code. 


! 

! 

I 
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